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Abstract. We consider the polarization properties of ra- 
diation emitted by relativistic charged particles while 
moving along the curved magnetic held lines in a pul- 
sar manetosphere. We propose that the radiation emitted 
by positrons and electrons while moving along the curved 
magnetic field lines is orthogonally polarized. The polar- 
ization angle of each orthogonal polarization mode is well 
described by the rotating vector model. However, the po- 
larization angle swings observed in micropulses and sub- 
pulses are often found to be contrary to the rotating vector 
model, and such swings are expected to arise due to the 
coherent superposition of orthogonal polarization modes. 

As an application of our model, we discuss the po- 
larization of pulsar PSR B0950+08 at the frequency 
1.71 GHz. Data on individual pulses, obtained using the 
100-m Effclsberg radiotelescope, was statistically analysed 
and the results are presented as probability of occur- 
rence gray-scale plots. We find that pulsars emit radiation 
mainly in the form of two independent orthogonal modes. 
It seems, they exist at all pulse longitudes but at some 
longitudes one mode dominates over the other. The polar- 
ization angle gray-plot indicates two most favoured angles 
separated by approximately 90° at each pulse longitude. 
The depolarization is mainly caused by the incoherent su- 
perposition of orthogonal modes. 

We infer from this study that pulsar radiation consists 
of two major elliptically polarized electromagnetic waves 
with orthogonal polarization angles. Our model predicts 
that such waves could be radiated by the positrons and 
the electrons accelerated along the curved magnetic field 
lines. The sense of circular polarization of these modes 
depend upon from which side of particle trajectory the 
radiation is received. 
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1. Introduction 



Pulsar radio emission and further propagation effects in 
the magnetosphere are not well understood. This situation 
is partly due to the difficulty in understanding particle 
acceleration and current flow in pulsars superstrong mag- 
netic field. From the observation point of view curvature 
radiation seems to be most attractive among all the pro- 
posed radio emission mechanisms (Michel 1991; Meszaros 
1992; Xilouris et al. 1994). However the coherent curva- 
ture emission by bunches has been often criticized (e.g. 
Melrose 1981). Curvature radiation is quite similar to syn- 
chrotron radiation, with the only difference that the role 
of Larmor radius is played here by the radius of curva- 
ture of magnetic field lines. Michel (1987a) has indicated 
that curvature radiation has significant circular polariza- 
tion when viewed at an angle to the plane of a magnetic 
field line, and reverses sense when viewed from the other 
side. 

The curvature radiation model has been developed by 
assuming that particles follow the curved field lines (Stur- 
rock 1971; Ruderman & Sutherland 1975). The theory of 
synchrotron radiation is developed by assuming uniform 
and straight field lines while the curvature radiation is 
developed by ignoring the spiral motion of the particles. 
The most often cited reason for ignoring spiral motion is 
that particles will instantly lose the energy associated with 
perpendicular component of motion by synchrotron radia- 
tion, cascading down to zero free energy. This is true near 
the polar cap where the field lines are straight. But when 
the particles move into the curved region of field lines 
they recover the perpendicular component of momentum 
as the field lines curve off from the direction of particle 
velocity. Hence the radiation emitted by charged particles 
while moving along the curved magnetic field lines cannot 
be described by considering either synchrotron or simple 
curvature radiation separately. Gil & Snakowski (1990) 
have attempted to examine the polarization properties of 
curvature radiation but they have also not considered the 
role of magnetic force. 
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One of the most fascinating features of pulsar radia- 
tion is the occurrence of orthogonal polarization modes 
(OPM), i.e., the two electromagnetic waves with orthog- 
onal electric vectors. This phenomenon has become diffi- 
cult to explain by emission models based on the simple 
curvature emission, as it cannot specify the two preferred 
polarization angles. There are other models based on the 
propagation effects but one would expect such effects to 
be strongly dependent on frequency and require special 
viewing angles (Michel 1991), whereas observations indi- 
cate the phenomenon is broad band. The interpretation 
based on the geometrical effects (Michel 1987b) indicates 
that two modes could be due to the overlap of radiation 
from two distinct emission regions in the magnctosphere. 
However, it is not clear why the two separated sources 
have nearly orthogonal polarization. 

Pulsars have been noted for their highly polarized ra- 
diation. To understand the phenomenology of pulsar po- 
larization, many attempts have been made to fit the av- 
erage polarization angle swing within the context of the 
rotating vector model (RVM) (Radhakrishnan & Cooke 
1969). However many pulsars do not fit with this inter- 
pretation (Manchester 1971; Rankin ct al. 1974) and such 
discrepancies have been attributed to the occurrence of 
OPM (Backer et al. 1976; Gil & Lyne 1995). The relative 
strengths of OPM vary with pulse longitude, thus caus- 
ing discontinuities in the average pulse polarization angle 
swing (Cordes & Hankins 1977). The switching between 
OPM mostly but not always occurs on the boundaries of 
micropulses and subpulses (Cordes & Hankins 1977), and 
may be intrinsic to the emission process (Manchester et al. 
1975; Cordes et al. 1978). The single-pulse studies (Stine- 
bring ct al. 1984a, b) show that OPM overlap, and in the 
vicinity of pulse longitude where jump occurs some pulses 
prefer to have one polarization and some have other polar- 
ization. All these fine details are lost when the pulses are 
averaged. The depolarization in average pulse occurs due 
to the superposition of OPM (Lyne et al. 1971; Manch- 
ester et al. 1975; Stinebring et al. 1984a). 

Cordes et al. (1978) have reported the existence of 
OPM with opposite senses of circular polarization in pul- 
sar PSR B2020+28 at 430 MHz. They found no evidence 
for a threshold intensity in the occurrence of OPM. The 
peak subpulse emission can be in either OPM, and some- 
times transitions between modes do not occur on the edges 
of subpulses. Therefore, the occurrence of OPM may be 
stochastic and perhaps indicating the untenability of a ge- 
ometric interpretation of the transitions as aspects of an 
angular beam of radiation. 

The motion of a charged particle along a rotating mag- 
netic field line is discussed by Gangadhara (1995; 1996a) 
and in the case of a curved magnetic field line (Gangad- 
hara 1996b), and the studies indicated that OPM may be 
produced by positrons and electrons, as they have got op- 
posite senses of gyration. In this paper, we show that the 
relativistic positrons and electrons can produce OPM be- 



cause while moving along a curved magnetic field line their 
accelerations become inclined with respect to the radius of 
curvature of a field line. We propose the coherent superpo- 
sition of OPM as a physical explanation for the polariza- 
tion angle swings observed in micropulses and subpulses. 
In Sect. 2 we discuss the forces acting on a positron and 
an electron while moving along a curved magnetic field 
line. The radiation fields due to a positron and an elec- 
tron are derived in Sect. 3. The polarization properties of 
radiation fields are computed in Sect. 4. Finally, in Sect. 5, 
we present a series of gray-plots describing the polariza- 
tion states of OPM from PSR B0950+08, and discuss the 
possible explanation based on our model. 

2. Motion of relativistic charged particles along 
curved magnetic field lines 

In the physics of pulsar magnctosphere, a very important 
question is how are the charged particles ejected from 
the surface of a neutron star, and it has been discussed 
by many authors (e.g. Herring & Nichols 1949; Good 
& Miiller 1956; Ruderman & Sutherland 1975; Kundt & 
Schaaf 1993). If particles are produced with the velocity 
having a component perpendicular to the magnetic field 
they quickly radiate away the energy associated with that 
component of motion by synchrotron radiation. Hence the 
particle gyration becomes almost absent near the polar 
cap. Once the particles radiate away the energy associ- 
ated with the perpendicular component of velocity, they 
are free to move along the field lines as long as they are 
straight. But when the particles enter the curved region 
of field lines, they recover perpendicular component of ve- 
locity at the expense of parallel component (Gangadhara 
1996b) as the field line curves off from the direction of 
particle motion (Fig. 1). This phenomenon must be true 
for both primary as well as secondary particles. 




B 



Fig. 1. Motion of a particle with charge q and velocity v along 
a curved field line B. 

The motion of a particle along a curved magnetic field 
line is governed mainly by the magnetic Lorentz force Fb 
and the centrifugal (inertial) force F c . Figure 2 shows the 
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directions of these forces when a positron and an electron 
are in motion along a curved magnetic field line. The net 
force acting on a particle is given by 



Fi = F ci + F 



Bi, 



where i = p for a positron and e for an electron, 
magnetic force is 



Bi 



Qi 

—v±i 

c 



x B 



and the centrifugal force is 



h Ci - -n P, 



(1) 



The 



(2) 



(3) 



where m\ and q\ are the mass and charge, 7; is the rela- 
tivistic Lorentz factor, v\\i and v±i are the components of 
velocity with respect to the magnetic field B, respectively. 
Here p is the radius of curvature of a magnetic field line. 




Fig. 2. The forces Fb p and F cp act on a positron, and Fee 
and Fee act on an electron while they are in motion along a 
curved field line B. Here p is the radius of curvature of field 
line. 



centripetal force to keep the particle on track. But the pe- 
culiarity of magnetic force is that it acts in the directions 
perpendicular to the magnetic field and not always in the 
direction of radius of curvature. Since \Fbi\ ~ \F c i\, the 
net forces F p and F c are inclined through, respectively, 
<~ 45° and —45° with respect to the plane of magnetic 
field line. Therefore, the net accelerations of positrons and 
electrons become orthogonal in the curved magnetic field 
lines of pulsar magnctosphere. The detailed computation 
of particle dynamics, including radiation reaction, is under 
consideration and published elsewhere. 



3. Radiation fields of positrons and electrons ac- 
celerated in curved magnetic field lines 

Consider a relativistic particle with charge q moving along 
a curved trajectory C (space curve) in the xyz-coordinate 
system (Fig. 3). Let 8 be the angle between particle po- 
sition r and yz-plane. A distant observer at P receives 
radiation at an angle 4> p from the plane of particle orbit. 
The electric field of radiation at the observation point is 
given by (Jackson 1976) 



fc x [(fc - g) x g] 
Sa 3 



(4) 



where a = 1 — /3 • k. The distance from radiating region 
to the observer is S, the propagation vector is k, and the 
velocity and acceleration of particle are (3 = v/c and (3. 

The radiation emitted by a relativistic charged particle 
has a broad spectrum. The range of frequency spectrum 
is estimated by taking Fourier transformation of electric 
field of radiation: 



+00 



E{u) = 



E{t)e iujt dt. 



(5) 



Initially, magnetic and centrifugal forces act in the di- 
rections perpendicular and parallel to the plane of mag- 
netic field line, respectively. The motion of particles along 
curved field lines is discussed by Jackson (1976) and Gan- 
gadhara (1996b). The gyration motion (spiral motion) is 
considered as the zero-order-motion, and the motion that 
arises due to the curvature of field lines as the first-order 
motion. In the absence of initial perpendicular component 
of velocity, particles cannot gyrate around the curved mag- 
netic field lines. The reason is that as particle attempts 
to spiral, field line curves off to the side and hence the 
centrifugal force prevents it from spiraling (Gangadhara 
1995; 1996a). Also, it is observationally known that no 
generalized Faraday rotation is evident in pulsar magne- 
tospheres (Cordcs 1983; Lyne & Smith 1990). There is no 
central force at O, and the magnetic force plays the role of 



In Eq. (4), ret means evaluated at the retarded time t' + 
^ — t. By changing the variable of integration from t 
to t 1 , we obtain 

E( U ) = 4= q - [ * X [{k ~ f } X ^ e M^(Q/c}< (6 ) 

—00 

where we have used dt = adt'. 

When the observation point is far away from the region 
of space where the acceleration occurs, the propagation 
vector k can be taken constant in time. Furthermore the 
distance S(t') can be approximated as 

S(t')nS -k.r(t'), (7) 
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where S is the distance between the origin O and the 
observation point P, and r(t') is the position of particle 
relative to O. Then Eq. (6) becomes 



loss of generality, to lie in a plane which is parallel to the 
xy-plane: 



k = (cos^p, sin0 p , 0), 



(12) 



+ 00 , 



EM « ^ / tx|(fe -/ )X ^^, (8) 



2-kSoC 



where the primes on the time variable have been omitted 
for brevity. The integrand in Eq. (8), excluding exponen- 
tial, is a perfect differential, therefore, we can integrate by 
parts, and obtain 



where P is the angle made by k with respect to the 
positron velocity. Using Eqs. (11) and (12), the vector part 
of the integrand in Eq. (9) can be written as 



k x (k x (3) = (3 



18 



cos < 



-^=-£|| +^ sin (/) p cos 9 + j i± 

(13) 



E(uj) = -i 



.quj e 



iuS /c 



-t-oo 



2^5*00 



k x (k x f3)c^ t - Lr / c Ut. (9) 



The polarization of emitted radiation can be estimated for 
a specified motion with known r(t) and (3. 



where i\\ = z is a unit vector in the direction of z-axis, 
and e± = — sin</> p x + cos (f> p y = i\\ x k is a unit vector 
which is orthogonal to both £\\ and k. 

Using Eqs. (10) and (12), the argument of exponential 
in Eq. (9) can be written as 



U)\t- 



k.r 



= U) 



R ( , . sin^pcos^ 

t cos <fin sin H -= 

c V V2 




(14) 



Since R is close to the radius of curvature of particle tra- 
jectory, 9 can be replaced by c(3t/R. Due to the relativistic 
beaming only for small values of <j) v will there be an ap- 
preciable radiation, and therefore the duration of pulse 
At <~ R/c-f becomes very short. Hence the arguments of 
sine and cosine functions in Eqs. (13) and (14) are of the 
order of I/7. So, we obtain 



fcx(fc x /3) w (3 



(3c f (3c 



1 



2i? 2 



Fig. 3. The geometry used to describe the radiation emitted 
by a positron while moving along a curved trajectory C. Here 
6 is the angle between yz-plane and r, and 4> p is the angle 
between (3 and k. 



and 



u)\t- 



k.r 



2 



(15) 



(16) 



Radiation field due to a positron 

Since the duration of pulse At ~ J?/c7 is very short, it is 
necessary to know the position r(i) and velocity (3 of par- 
ticle over only small arc of the trajectory whose tangents 
are in the general direction of observation. Therefore, for 
a positron moving under the action of magnetic and cen- 
trifugal forces, we take: 



r(t) — R (sint 



cos 9 cos 9 



V2 ' V2 



(10) 



where (3 has been set equal to unity at wherever possible, 
and neglected all those terms which are of the order of 
I/7 2 times those kept. 

The components of electric field (Eq. 9) in the direc- 
tion of unit vectors eii and are: 



Jt exp 



1 



2 I V 7 

2 



t + 



dt, (17) 



/3(f) = (3 (cos9, 



sin 9 sin 9 



V2 ' V2 



(11) E ± (w) 



Since the integral in Eq. (9) has to be computed over the 
path of particle, the unit vector k can be chosen, without 



c 2 (3 2 
2i? 2 ' 



So/C [ ljL t . 

2nS c J I V2R 



-1 , e 



exp 



.U) 



1 



+ <t>; )t + 



3i? 2 



dt. (18) 
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in sin 9 N 



To reduce the integrals into some known forms, we change 
the variable to 



(3{t) = (3 cos 



(30) 



y = 



<P 2 P 



-1/2 



Ct 

R 



and introduce a parameter 



ujR f 1 



3/2 



3c \-f z 
Therefore, the integrals are 



(19) 



(20) 



V2 ' V2 

We choose the propagation vector for electron field as 

k = (cos C , sin C , 0) (31) 

such that it is parallel to the k in Eq. (12) and ob- 
server receives the radiation from both particles (positron- 
electron) . Here e is the angle between k and the electron 
velocity. Using Eqs. (30) and (31), the vector part of the 
integrand in Eq. (9) can be written as 



1 



' -(- 



1/2 



+ 00 



— \\dy, (21) kx{kxf3) = (3 



h = 



h = 



R 2 



R :) 



2+€ 



— oo 



smfc 

7f 



■en + sin O cosd - 



cos 0o sin ( 



V2 



V exp<j i-£ p ( y 



3 



| dy, (22) wnere — z and e±_ = — sin0 c x + cos0 o y = £|| x k. 

The argument of the exponential in Eq. (9) is 



(32) 



2 +K 



3/2 



J y 2 CXp|i; 



y + 



dy. 
(23) 



sin O cos 
V2 



(33) 



We can identify these integrals with Airy integrals, and Using the identical approximations and solutions of inte- 
grals used in the previous section for a positron field, we 
1/2 find the components of an electron field: 



obtain the solutions 
2 R ( 1 



h = 
h = i 



V3 c V7 2 
2 R 2 ( 1 

2 



V3 cr 



7 



Ki/s^p), 
K 2 / 3 (Cp), 



Ja = - 



2 i? 3 / 1 



y/Z c 3 V7 



v2 + n 



3/2 



Kl /3 «p), 



(24) 
(25) 
(26) 



1 q(i 2 UjR SJ 
^||c(w) = — ^= — rr^e 



where K4/3 and K 2 /3 are modified Bessel functions. Sub- 
stituting the solutions of integrals into Eqs. (17) and (18), 
we get 

qf3 2 tuR 



3tt c 2 ^ 
YqPuiR 



' (^+0 e 2 ) K 2 :i (< ( .). (34) 



3tt c 2 S q 
1/2 



1 



^5 /c 



/? 2 



V2V7 2 



K 2/3 (Cc) 
Ki/ 3 (Co) 



£|| P M = 



3irc 2 S 



1+0 2 )k 2/3 (£ p ), (27) 



(35) 



4. Polarization of radiation emitted by positrons 
and electrons 



E ±p (uj) 



2 qfiuR ^sjc 

3-7T C 2 ^ 



JL (1_ 

V2W 



— + ' 

1 2 



1/2 



The electric fields E p (uj) and E c (uj) derived in Sect. 3 de- 
scribe the polarization properties of OPM. It is clear that 
one mode is emitted by positrons and other by electrons. 
(j) 2 ) )■ K 1 / 3 (^ p ) , If the two radiation fields do not bear any phase relation 
then they are expected to be incoherently superposed at 
(28) the observation point. On the other hand, if there is a 

phase relation then they are coherently superposed. From 
where we have introduced a suffix 'p' on E\ { and E ± to thfi observational point of view both the cases are impor _ 



K 2/3 (e P 
Ki/3(£p) 



indicate that they are due to a positron. 

3.2. Radiation field due to an electron 

For an electron, the position and velocity can be taken as 



r(t) =R[sm0,- 



cos cos 6 



V2' V2 



(29) 



tant, and we discuss them separately in the following two 
subsections. 

4-1. Incoherent superposition of radiation fields 

Consider a region in the magnetosphere containing a large 
number of radiating positrons and electrons. Let each par- 
ticle emit a pulse of radiation with electric field E (t). An 
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observer will detect a series of such pulses, all with same 

shape but random arrival times t\, £2, ts, £n- Then the 

measured electric field will be (Rybicki & Lightman 1979) 



N 



E(t) = J2 E o(t-h), 

j=l 



where N is the number of particles. Taking the Fourier 
transform, we find 



1 " r-t-00 
Let u = t — t j then du = dt. Therefore, we get 



(37) 



N 



E[u) = £ (w)^V w *J. 

j=i 



(38) 



Let dW, dA and du be the differential incriments in 
energy, area and frequency, respectively. Then the mea- 
sured spectrum is given by 



dW 
dAduo 



c\E(co)\ 2 

= C |£ Mi 2 f;fy^-^ 
j=i k=i 



= c\E {uj) 



N 



N + cosw(tj — <k) 



Since tj and ik are randomly distributed in the case of 
radiation fields which do not have any phase relations, 
the second term averages to zero. Therefore, we have 



dW 
dAduo 



c\E(lu)\ 2 = Nc\E q (w)\ 2 . 



(40) 



Hence in an incoherent superposition of radiation fields 
the measured intensity will be simply a sum of intensities 
radiated by the individual charges. 

4-2. Coherent superposition of radiation fields 

The very high brightness temperature (10 25 -10 30 K) of 
pulsars, lead to the conclusion that the radiation must 
be coherent. Pacini & Rees (1970), and Sturrock (1971) 
among others were quick to point out that the observed 
coherence may be due to bunching of particles in the 
emission region of magnetosphere. However, the topic of 
bunching mechanism continues to be an outstanding chal- 
lenge. If the bunches of plasma particles with sizes much 
smaller than a wavelength exist then the arrival times 
tj i=a 0, because all the pulses will have a same arrival time 
to order (size of bunch) / (wavelength) . Then Eq. (38) gives 



E(w) « NE (lu). 



(41) 



Hence the total radiation field due to a bunch of particles 
is equal to the vector sum of the fields radiated by each 
charge. 

Now, the measured spectrum is given by 



(36) dW 



dAduj 



= c\E(uj)\ 2 =N 2 c\E (lu) 



(42) 



Hence the coherent sum of radiation fields of a bunch of 
particles leads to the intensity, which is equal to iV 2 times 
the intensities due to the individual charges. 

Since the pair creation, breakdown of polar gaps and 
sparks are not steady state processes, the plasma in the 
bunches can be neutral or nonneutral at an arbitrary time. 
If the bunches consist of prominently positrons (electrons) 
then the radiation field will be a coherent sum of E p of 
each positron (E c of each electron). In such cases, the 
two radiation fields are orthogonally polarized, and the 
modes are said to exist disjointly. On the other hand, 
if bunchs are neutral then both charges can contribute 
equally and simultaneously to the radiation field. Since 
both the modes are coherently superposed in this case, 
we call this case joint. On the basis of individual pulse 
polarization Stinebring et al. (1984a) have concluded that 
OPM are superposed. If OPM exist disjointly their po- 
larization properties are well defined, but not when they 
exist jointly. In the following two subsections we consider 
these cases in detail. 



(39) 

4.3. Presence of OPM disjointly 



The OPM are said to exist disjointly when only one mode 
exists or one mode has become very strong compared to 
other mode. In all such cases the polarization state is solely 
determined by the stronger mode. The polarization state 
of OPM can be described more accurately using the Stokes 
parameters: 

h = EuETu + E ±i E*l ; 



{ 



1)1 

1 



2 



K5/s(6) + #(^ 



x<^ 1 + 



Qi — E\\iE*^ - E±iEj_i 



K 2 /3 te) 



2 



(43) 



Ui = 2Re[K i £ ±i ] 



\7 2 



K 2 /3 te), 
(44) 



(45) 
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Vi = 21m[E^E ±i ] 

= -V2a 2 fi<i 



3/2 



xK 1/3 te)K 2/3 te), 



where 



ujR ( 1 

Ik w 



6 = — 



3/2 



/3 2 / 1 



2 



(46) 



3tt c 2 5 g 



i = p and r\ — +1 for positrons, and i = e and 77 = —1 
for electrons. The intensity and polarization angle of the 
linearly polarized radiation are: 



L i = y fuf+Q*, 



1 -1 fUi\ 



(47) 



(48) 



Due to the action of forces Fbi and -F c i, particle orbital 
planes become inclined with respect to the plane of mag- 
netic field line. The radiation beams of the two charges 
appear on either sides of the plane of magnetic field line. 
Therefore, we define <f> v = (f> — - and <j> e = <f>+ -, where 4> 

is the angle between the plane of magnetic field line and k. 
Using 7 = 300, v = 1 GHz and R = 10 6 cm, we computed 
the polarization parameters h, Li, Vi and tpi, and plotted 
as functions of <f> in Fig. 4. The continuous line curves de- 
scribe the positron radiation field while the broken line 
curves for an electron radiation field. About the particle 
orbital planes, i.e., about 4> p = and <f> e = 0, intensity and 
linear polarization have maxima while circular polariza- 
tion undergoes sense reversal. Therefore, when observed 
with line-of-sight lying in the range — — < <j> < — as indi- 
cated by a two way arrow in Fig. 4(c), one mode becomes 
right hand circularly polarized (V p > 0) while the other 
becomes left hand circularly polarized (V c < 0). Individ- 
ual pulses from PSR B2020+28 (Cordes et al. 1978) and 
PSR B0950+08 (Fig. lOg) indicate OPM tend to have op- 
posite circular polarization. Figure 4(d) shows the polar- 
ization angle as a function of <fi, clearly, two modes are or- 
thogonally polarized when Ii and Li are in maxima. Since 
particles are constrained to follow the curved field lines, 
the polarization angle of each mode swings in accordance 
with the RVM. In the case of PSR B0329+54 at 408 MHz, 
Gil & Lync (1995) have clearly shown that each of these 
modes is well described by the RVM. 

The Stokes parameters of E v {uo) and E c (u) clearly in- 
dicate OPM are elliptically polarized. The radiation field 



E p is polarized with ~ 45° while E c with 



-45° with 



respect to the plane of magnetic field line. Figure 5 illus- 
trates the polarization ellipses of OPM in the coordinate 
system with axes along the unit vectors fy, e± and k, 
where the unit vector is parallel to the radius of curva- 
ture (p) of magnetic field line. 
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Fig. 4. Polarization parameters versus <f> for the radiation due 
to a positron (continuous line curves) and an electron (bro- 
ken line curves). The labels on vertical axes of panels (a), 
(b) and (c) are dimensionless as we have normalized with 
p? = [g/3 2 a.7?/(3^) 1 / 2 c 2 S' 7 2 ] 2 . 



As mentioned in Sect. 4.2, particle generation pro- 
cesses such as direct emission of electrons from the stellar 
surface due to the building up of high potential differ- 
ence at the polar cap and the pair creation are not steady 
processes. They operate on time scales less than microsec- 
onds, therefore, the plasma which flows along the field 
lines is not always neutral. On time scales such as sam- 
pling time, it may look either negative or positive. In all 
such cases, the observer receives OPM at ~ ±45°, depend- 
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ing upon the sign of plasma which exist at that particular 
instant. One mode is not radiated even when the plasma 
is neutral, if the energy of one kind of charge happens to 
fall below the threshold to radiate. 





V 45 ° 






A 

k 





Fig. 5. Representation of OPM in a coordinate system with 
axes parallel to the orthogonal unit vectors eii , ej_ and k, where 
E|| is parallel to the radius of curvature (p) of magnetic field 
line. 



4-4- Presence of OPM jointly 

If the charges are close enough to each other (e.g. bunch), 
we may expect the fluctuations in amplitudes and phases 
of radiation fields will not be independent at the obser- 
vation point. In the limit where the distance between the 
two charges is much less than a wavelength, the ampli- 
tudes vary in phase (Born & Wolf 1986). Therefore, when 
the Fourier components E p (oj) and E c (u>) are coherently 
superposed, we get 

E(u) = E p {u) + E e (w). (49) 

The Stokes parameters of total field can be defined sim- 
ilar to those (Eqs. 43-46) for OPM. Again using 7 = 300, 
v = 1 GHz and R = 10 6 cm, we computed polarization 
parameters I, L, V and ip of E(w), and plotted as func- 
tions of 4> in Fig- 6. The two humps at <~ 0.2° and —0.2° 
in Figs. 6(a) and (b) are due to the dominance of emission 
from positrons and electrons at those angles, respectively. 
Figure 6(c) shows the sense reversal of circular polariza- 
tion about — 0: i- e -> about the plane of magnetic field 
line. The polarization angle of E{lo) is plotted in Fig. 6(d), 
and the swing arises due to the coherent superposition of 



OPM. This type of polarization angle swings are quite 
evident in micropulses and subpulses. 

Consider a coordinate system- xyz such that the major 
axes of polarization ellipses of E p (t) and E c {t) are parallel 
to the axes x and y, respectively, as shown in Fig. 7. For 
the sake of illustration of OPM we represent them as 

E p (t) = E xp cos(u>i + 5 p )x + E yp sm(ujt + 5 p )y, (50) 

E c (t) = E xc sin(wt + 5 e )x + E yc cos(ui + S e )y, (51) 

where _E X ;, E y \ and <5; are the amplitudes and initial 
phases. The period of rotation of electric fields is T = \ jv. 
For v = 1 GHz we find T = 10~ 9 s, much smaller than the 
sampling interval in a typical observation (> 1/xs). 

Let S = dp — S e be the phase difference between OPM, 
and tp be the polarization angle of total field E(t) mea- 
sured counter-clockwise from the x-axis. In Fig. 7, the el- 
lipse drawn with dots indicate the polarization state of 
E(t) with the sense of rotation can be clockwise or coun- 
terclockwise or even nil, depending upon strengths of cir- 
cular polarization of OPM. For the values of 5 in the range 
— I to ^ we find ip lies between and |. On the other 
hand for values of S in the range ^ to 4p , ip takes values 
between | and it. The polarization state of E{t) follows 
closely the stronger mode. 

It is probable that separate streams of positrons and 
electrons exist in the pulsar magnetosphere (Sturrock, 
1971). The mode changing observed in PSR B0329+54 
by Hesse ct al. (1973) indicate the channeling of particles 
along field lines is not steady. Assume that there exist two 
separate and closely spaced streams: one consists of mainly 
positrons and other electrons at longitudes 0i and 02, re- 
spectively. When the observer's line-of-sight moves from 
0i to <p2 there will be an orthogonal jump in the polariza- 
tion angle. Now, what happens to degree of polarization 
during those jumps? Let St be the temporal separation be- 
tween 0i and </>2- When one approaches 2 starting from 
0i, field E p becomes progressively weaker while E e gets 
stronger. To find the polarization angle distribution during 
St, let us divide the interval St into n subintervals having 
the widths equal to the period (T = \jv) of rotation of 
E(t). If 5 lies between — | and | then E(t) will have el- 
lipses t\ =T, t 2 = 2T, .... t n = nT = St with polarization 
angles ip2, ■■■■ ipn, respectively, as shown in Fig. 8. On 
the other hand when S lies between ^ and ^ we get the 
ellipses t\ = T, r 2 = 2T, .... t„ = nT = St with polar- 
ization angles —tpi, —1P2, ■■■■ —tpn- If the observations are 
made with the time resolution, which is of the order of St, 
we get total depolarization. This is because during the in- 
tegration the electric fields oriented in all directions (doted 
line ellipses) are superposed. The general result of obser- 
vation that the orthogonal polarization angle jumps are 
accompanied by percentage of polarization going to zero, 
are caused by such an effect. For example, see the pulses 
of PSR B1804-08 and PSR B1905+39 by Xilouris et al. 
(1991). 
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Fig. 7. Representation of OPM and total electric field E(t) in 
the coordinate system-xyz. The angle ip represent the polar- 
ization angle of total field. 
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Fig. 6. Polarization parameters versus <j> for the total radi- 
ation field E(w). The labels on vertical axes of panels (a), 
(b) and (c) are dimensionless as we have normalized with 
= [qp 2 cjR/(3^ 2 c 2 S o7 2 ] 2 . 



5. Discussion with an application to polarization 
of PSR B0950+08 

In the previous sections we discussed the polarization 
properties of radiation emitted by positrons and electrons 
while moving along a curved magnetic field line. Here we 
extend our discussion to the emission from plasma stream- 
ing along many field lines. In the curved region of magnetic 
field lines, for the reasons given in Sect. 2, positrons move 




Fig. 8. Orthogonal polarization angle jumps occur when one 
mode say E p becomes weaker while the other mode E c gets 
stronger or vice versa over a small interval of pulse phase during 
pulsar rotation. 
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to one side of field line while electrons to other side as 
shown in Fig. 9. Let us consider three observing positions: 
A at the left side of the pulse, B at the middle and C at 
the right side. 
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Fig. 9. Projected field lines loaded with positrons (+ signs) 
and electrons (— signs). Points A, B and C are the three ob- 
serving positions. 



At the position A, observer mainly receives radiation 
from the left side of planes of magnetic field lines, and the 
circular polarization of OPM tend show opposite senses 
(see Fig. 4c). An observer at position B (fiducial plane) 
views particle trajectory planes edge-on, therefore, the cir- 
cular polarization of each mode undergoes sense reversal. 
Finally, an observer at C receives radiation mainly from 
the right side of planes of magnetic field lines, and circular 
polarization of each mode becomes opposite to its value 
in position A. If the OPM exist disjointly then the polar- 
ization angle of each mode swings in accordance with the 
RVM when the observers line-of-sight moves from point 
A to C. 

Radhakrishnan & Rankin (1990) have identified the 
two extreme types of circular polarization signatures: an 
antisymmetric type (Fig. 6d) wherein the circular polar- 
ization changes sense at the middle of the pulse, and a 
symmetric type wherein it is predominantly of one sense. 
Our model can easily explain the antisymmetric type of 
circular polarization, while for symmetric type we need 
further analysis on the high resolution data of individual 
pulses. 



5.1. Observations of PSR B0950+08 

The pulsar PSR B0950+08 was observed in 1994 April 
by using the 100-m Effelsberg radiotclescope. Using a 
tunable HEMT-receiver with a system temperature of 
28 K, observations were made at the center frequency of 
1.71 GHz with a bandwidth of 40 MHz. The two circu- 
lar polarizations are separated in the receiver and am- 
plified. The signal is then fed into an adding polarime- 
ter, a passive device with four output channels which al- 
lows further online signal processing. The pulse-smearing 
caused by the dispersion due to the interstellar medium 
is than removed using an online dedispersion device. This 
is a four unit 60 x 667 kHz filterbank. The output of each 
channel is then detected and converted into a digital sig- 
nal by a fast A/D converter. After a time delay accord- 
ing to the dispersion measure, the outputs of all chan- 
nels are added and than recorded by the backend. Af- 
ter a careful calibration procedure, Stokes parameters are 
obtained from the four recorded output channels. A de- 
tailed system description and the calibration procedure 
are given by von Hoensbroech & Xilouris (1997). At fre- 
quency 1.71 GHz, we recorded about 1200 pulses with a 
time resolution of 0.24 ms. 

We consider only the main pulse as the interpulse is 
too weak to reproduce the properties of OPM. The average 
polarization parameters are plotted as functions of pulse 
phase <p in Fig. 10(a). The continuous line curve indicates 
intensity (7) variation while broken and doted ones indi- 
cate linear (L) and circular (V) polarization, respectively, 
in arbitrary units. 

The gray-plots show the frequency of occurrence of 
OPM at different pulse phases, and have become power- 
ful tools in analyzing the pulsar polarization properties. 
The darkest regions represent the most probable regions 
of occurrence. For each pulse phase bin where I is above 
3tr level of its value in the off pulse region, gray-plots were 
computed for L, V, I and polarization angle ip. All those 
phase bins, where the condition L 2 + V 2 < I 2 was not met, 
were excluded as they lead to spurious interpretation of 
polarization. The points with error bars in Fig. 10(b) indi- 
cate the integrated polarization angle superposed over the 
polarization angle gray-plot. Figures 10(c)-(e) represent 
the gray-scale maps of percentage of linear polarization 
[L(%) = 100L/7], circular polarization [V(%) = 100V/I] 
and I, respectively, while Figs. 10(f)-(h) represent the 
gray-scale scatter plots of L(%), V(%) and I versus the 
polarization angle ip. 

The polarization angle gray-plot (Fig. 10b) shows the 
two most preferred tracks close to 45° and —45° posi- 
tion angles in the ip versus <p plane, in agreement with 
our model (Fig. 4d). The frequency of occurrence of OPM 
with respect to the polarization angle is shown in Fig. 11, 
proving the importance of OPM in the pulsar radiation. At 
any pulse longitude the average polarization angle curve 
follows closely the mode which is more intense, as indi- 
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Fig. 10. Polarization histograms of pulsar PSR B0950+08 at frequency 1.71 GHz: (a) average pulse with arbitrary units (a.u), 
(b) polarization angle if), (c) & (f) percentage of linear polarization, (d) & (g) percentage of circular polarization, and (e) & (h) 
intensity. The angle <f> denotes the pulse phase. 
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cated by Fig. 10(b). The uniform or random component 
of polarization angle is probably due to the coherent su- 
perposition of OPM (see Figs. 6d&7). This idea is also 
supported by the observation that the random compo- 
nent becomes significant only at those pulse longitudes 
where both modes exit. For example, see the polariza- 
tion angle displays of PSR B0823+26, B0950+08 (Stinc- 
bring et al. 1984a), B0834+06 (Stincbring et al. 1984b), 
and B0329+54 (Gil & Lync 1995). McKinnon and Stinc- 
bring (1996) have also suggested that the random compo- 
nent may arise from superposed modes. Barnard & Arons 
(1986) have proposed lo~ 2 frequency dependence for an- 
gular separation of OPM, but our model do not predict 
any appreciable frequency dependence. 
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Fig. 11. Frequency of occurrence of pulses with respect to the 
polarization angle. 
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On the basis of relative frequency of occurrence 
of OPM at different pulse longitudes as indicated by 
Figs. 10(b), we may identify pulse longitude ranges I, II 
and III, as marked in Fig. 10(a). It is clear that both 
modes exist at all pulse longitudes but one mode dom- 
inates over the other at some longitudes. For example, 
in region II the mode (say mode 1) with tp <~ 45° dom- 
inates while in region III the other mode (mode 2) with 
ip ~ —45° dominates. 

When the frequency of occurrence of OPM is very high 
as in region II, it is likely that more often they simulta- 
neously exist. If so then they get superposed either co- 
herently or incoherently depending upon whether there 
is any phase relation between them or not. Figure 10(d) 



Fig. 12. Grey-scale plots representing the correlations be- 
tween L(%), V(%) and I in the individual pulses. 



shows, at any pulse longitude, circular polarization can be 
of either sense depending upon from which side of parti- 
cle trajectory the radiation is received (see Fig. 4c). The 
distribution of different pulse intensities (arbitrary units) 
versus 4> is shown in Fig. 10(c). 

In revealing the properties of OPM, Figs. 10(f)-(h) 
are much more expressive than Figs. 10(c)-(d) as the ear- 
lier ones directly display the distributions against the po- 
larization angle if). Figures 10(f) and (g) indicate that if 
the modes exist disjointly, they appear at <~ ±45° with 
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high inear and circular polarization, but the random com- 
ponent is less polarized. Figure 10(g) shows OPM tend 
have opposite circular polarization (see Fig. 4c&d). Fig- 
ure 10(h) shows the distribution of intensities of OPM 
versus the polarization angle, and it indicates the mode 1 
is stronger than the mode 2. 

From the theoretical point of view, it is very impor- 
tant to find the correlations between L{%), V(%) and I. 
Using the individual pulse data, we computed the grey- 
plots (Fig. 12), which represent the correlations between 
different polarization parameters. Figure 12(a) shows lin- 
ear polarization becomes maximum when the circular po- 
larization is at minimum, a prediction of curvature radi- 
ation (see Figs. 4b&c). 

The behaviour of L(%) with respect to / is shown in 
Fig. 12(b). It shows anticorrelation between / and L(%) 
at higher intensities. The sharp cutoffs close to the verti- 
cal axis is due to the condition that / is above 3cr level. 
Manchester et al. (1975) and Xilouris et al. (1994) have 
also indicated the anticorrelation between / and L{%). 
The behaviour of circular polarization with respect to the 
intensity is shown in Fig. 12(c). At higher intensities, cir- 
cular polarization is lower. The reason could be the super- 
position of OPM with opposite senses of circular polariza- 
tion. 

6. Conclusion 

The motion of particles along the curved field lines can- 
not be described in analogy with the motion of particles 
in a central force. The radiation emitted by positrons and 
electrons while moving along the curved magnetic field 
lines is probably orthogonally polarized. The radiation is 
highly polarized when the OPM are not superposed. The 
polarization angle of each mode swings in accordance with 
the RVM. However, the polarization angle swings observed 
in micropulses and subpulses are often found to be con- 
trary to the predictions of RVM. We expect such swings 
are produced when the OPM are coherently superposed. 
Our model do not predict any appreciable frequency de- 
pendence on the angular separation between OPM. The 
circular polarization of OPM tend to have opposite senses. 

The observations of PSR B0950+08, particularly, po- 
larization angle, linear and circular polarizations of OPM 
are in agreement with our model. The polarization his- 
tograms clearly indicate that the depolarization is due to 
the superposition of OPM. This effect becomes much more 
severe at higher frequencies as the coherence-length which 
is of the order of wavelength becomes small leading to an 
incoherent superposition of OPM. 
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